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ABSTRACT 



Context. The majority of observed mass-to-light ratios of globular clusters are too low to be explained by 'canonical' 
cluster models, in which dynamical effects are not accounted for. Moreover, these models do not reproduce a recently 
reported trend of increasing M/L with cluster mass, but instead predict mass-to-light ratios that are independent of 
cluster mass for a fixed age and metallicity. 

Aims. This study aims to explain the M/L of globular clusters in four galaxies by including stellar evolution, stellar 
remnants, and the preferential loss of low-mass stars due to energy equipartition. 

Methods. Analytical cluster models are applied that account for stellar evolution and dynamical cluster dissolution to 
samples of globular clusters in Cen A, the Milky Way, M31 and the LMC. The models include stellar remnants and 
cover metallicities in the range Z = 0.0004 — 0.05. 

Results. Both the low observed mass-to-light ratios and the trend of increasing M/L with cluster mass can be reproduced 
by including the preferential loss of low-mass stars, assuming reasonable values for the dissolution timescale. This leads 
to a mass-dependent M/L evolution and increases the explained percentage of the observations from 39% to 92%. 
Conclusions. This study shows that the hitherto unexplained discrepancy between observations and models of the mass- 
to-light ratios of globular clusters can be explained by dynamical effects, provided that the globular clusters exhibiting 
low M/L have dissolution timescales within the ranges assumed in this Letter. Furthermore, it substantiates that M/L 
cannot be assumed to be constant with mass at fixed age and metallicity. 
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1. Introduction 



The mass-to-light 
have been given 



ratios 
a lot 



of 
of 



globular clusters (GCs) 
atte ntion recently (e.g. , 
McLau ghlin fc van der Marell '2005: Reikuba et"all l2007t 



Mieske fc Kroupal 120081 : Il5abringhausen et all 120081 ) 



Reikuba et al. ( 20071 ) have observed an M/L trend with 



cluster mass above a certain cluster mass, as more mas- 
sive clusters appea r to have higher M/ L than low-mass 
clusters (see also iMandushev et al.l Il99l[ ). This is an 
observation contrary to fundamental plane studies of 
GCs (e.g., iMcLauehlinl [200l and also in strong dis- 
agreement with the constant M/L for fixed age that 
is commonly assumed in observational a nd theoret- 
ical GC studies (e.g., [Harris et aT] 120061 : iMora et al] 
I2007t iBekki et all 120071) . Moreover, for Galactic GCs 
IMcLauehlinl (|2000l ) find M/Ly = 1.45 Mp, L;;\ whereas 
Simp l e Stellar Population models (e.g ., Bruzuai fc Charlo"3 



I2003t lAnders fc Fritze-v. Alvenslebenil2003D predict signif- 
icantly higher values of M/Ly = 2 — 4 Mq L^ for typical 
GC metallicities. Given the important role of GCs in 
galactic astronomy, it is essential to explain these apparent 
contradictions. 

In numerical and analytical stu dies of dynamical effects 
in clu s ters (e.g.,[B aumgardt fc M akinoll2003l : iLamers et all 
[200l IkT^ ssen fc LamersI 120081) it has become clear 
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that the dynamical evolution of clusters strongly af- 
fects cluster luminosity, co lour and mass-to-light ratio. In 
iKruijssen fc LamersI (|2008l hereafter KL08) it is shown how 
the evolution of these observables changes due to dynamical 
effects such as the preferential loss of low-mass stars and 
the retain of stellar remnants, but also due to the stellar 
initial mass function and metallicity. It is shown that M/L 
cannot be assumed to be constant for a fixed age and metal- 
licity, but instead depends on cluster mass when dynamical 
effects are accounted for. 

In this Letter, the analytical cluster models from KL08 
are applied to ex plain the obse r vation s of GCs in sev- 
eral g alaxies from iReikuba et all (|2007l ) and iMieske et al.l 
(|2008f ). In Sect. |5|I first summarise the models presented 
in KL08, which is applied to the observations in Sect. [3] In 
Sect. 13. li the effect of metallicity and the cluster dissolution 
timescale on cluster evolution in the {M, M/Ly}-plane is 
investigated. The observations are discussed in Sect. 13.21 
and are compared to the models in Sect. 13.31 A discussion 
of the results and the conclusions are presented in Sect.jH 



2. Modeling method 

In this study, analytic cluster models are used that incorpo- 
rate the effects of stellar evolution, stellar remnant produc- 
tion, cluster dissolution and energy equipartition. They are 
summarised here and are treated in more detail in KL08. 
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In the models, clusters gradually lose mass due to stellar 
evolution and dissolution. The total cluster mass evolution 
is described by dMd/dt — (dAfci/dt)cv + (dMci/di)dis, with 
the first term denoting stellar evolution and the second rep- 
resenting dissolution. Taking into account the formation of 
stellar remnants and the mass-dependent loss of stars by 
dissolution, this provides a description of the changing mass 
function and cluster mass in remnants. 

Stellar evolution is included by using the Padova 1999 
isochrone^. It removes the most massive stars from the 
cluster and increases the dark cluster mass by turning stars 
into remnants, which is i ncluded by as suming an initial- 
remnant mass relation. A iKroupal (|200l[ ) IMF is assumed. 

Cluster dissolution represents the dynamical cluster 
mass loss due to stars passing the tidal radius. This mass 
loss acts on a dissolution timescale Tdig: 



dMeift) 

dt 



Tdis 



to 



(1) 



where Tdis is related to the present day cluste r mass 
Mci(t) as Tdis = toMciity (|Lamers et ai.l 12005a'). lead- 
ing to the second equality in Eq. [TJ The characteristic 
timescale to depends on the environment and determines 
the strength of dissolution. For example, in the case of 
tidal dissolution to depends on tidal field strength and 
thus on galactocentric radius and galaxy mass. Typical 
values are = 10^ — 10* yr (e.g.. iLamers et aP [2005b) . 
translating into a total disruption time id°s'*' ~ 10* — 
10^"'^ yr for a 10^ Mq cluster. From TV-body simulations 
of tidal dissolution (ILa mers e t al. 2005^ and o bservations 
(jBoutloukos fc LamCTdl2003tlGieles et al.ll2005[ ). the expo- 
nent 7 is found to be 7 « 0.62. 

The effect of dissolution on the mass function depends 
on the dynamical state of the cluster. If it has reached 
energy equipartition, i.e., after core collapse, the cluster 
becomes mass-segregated and low-mass stars are preferen- 
tially los t. This occurs at about 20% of the total cluster 
lifetime (jBaumgardt fc Makinol [2003f) . For clusters with- 
out equipartition, bodies of all masses are lost with similar 
probabilities. 

Cluster photometry is computed by integrating the stel- 
lar mass function over the stellar isochrones, yielding clus- 
ter magnitude evolution M\{t, Md^i) for a passband A and 
a cluster with initial mass Md i. 



3. Applying the cluster models to observed clusters 

In this section I present the evolution of clusters in 
the {M, M/Ly }-plane, and apply this to explain the 
{M, M/Lv} distribution observed in real clusters. 



3.1. Cluster evolution in the mass-M/Lv plane 

In 'canonical' Simple Stellar Population (SSP) models, clus- 
ters only evolve due to stellar evolution, and therefore their 
mass-to-light ratios do not change due to dynamical effects. 
As the most massive stars (with low M/Ly) gradually dis- 
appear, in these models M/Ly is a monotonously increas- 
ing function of time that is constant for any set of clusters 
at a single age and metallicity. However, this is only correct 



if cluster dissolution occurs independently of stellar mass 
and the shape of the stellar mass function is preserved, i.e., 
there is no preferential loss of low-mass stars. 

Contrary to clusters from SSP models, in re- 
ality clusters do preferentially lose low-mass stars 
fe.g.. iHillenbrand fc HartmannI [l998t lAlbrow et all 120021 : 
iBaumgardt fc Makinoll2003f l. KL08 have shown that this 
strongly affects the M/L evolution of clusters due to the 
preferential loss of low- mass stars (having high M/L), and 
that consequently M/L cannot be assumed to be constant 
for clusters of a given age. 

In order to explain the appearance of clusters in the 
{M, M/Ly}-plane, I first compute the model cluster evo- 
lutionary tracks for different metallicities and dissolution 
timescales. The results are shown in Fig. [Tj 

Initially, model clusters are not in equipartition and 
they evolve to lower masses and increasing M/Ly (due 
to the death of massive stars), moving up and to the left 
in Fig. [TJ From the moment they reach mass segregation, 
happenin g earlier for low-mass cluste rs due to quicker re- 
laxation (jBaumgardt fc Makinol 120031 ) . they preferentially 
lose low-mass stars, which have high M/L. This decreases 
the cluster M/L and explains the maximum in the cluster 
evolution curves for lower initial masses. Along the cluster 
isochrones of constant age (dotted lines in Fig.[T|), M/Ly in- 
creases with mass, since at any given age low-mass clusters 
have spent more time in energy equipartition and thus have 
retained more massive (i.e., low-M/L) stars compared to 
massive clusters. The cluster isochrones flatten at the high- 
est masses, since these clusters have yet to reach equiparti- 
tion, leaving them at constant M/L. 

Within a galaxy, its GC s generally have similar ages 
fe.g.. lVandenberg et a l.' 1990*). Observations of GC systems 
should thus approximately follow cluster isochrones in the 
{M, M/Lv}-plane. Therefore, the isochrones are used in 
Fig. [1] to probe the influence of metallicity and dissolu- 
tion timescale on the expected GC distribution. Increasing 
metallicity leads to a higher maximum M/L, and thus also 
to steeper isochrones. Increasing the dissolution timescale 
shifts the entire isochrone to the left: for long dissolution 
timescales, it takes more time to reach equipartition and 
therefore only clusters of the lowest masses have preferen- 
tially lost low-mass stars. The dissolution timescale thus 
sets the location of the 'knee' in the isochrones, which is 
the cluster mass at which they flatten due to the absence 
of equipartition. This can be used to determine the disso- 
lution timescale range of an observed GC system. 



3.2. Observations of globular cluster mass-to- light ratios 



In iReikuba et~al] (|2007l . Fig. 9) the (dynamical-) 
{M, M/Ly}-plane is presented for GCs in several host 
galaxies. In Fig. [21 I show their {M, M/Lyj-plane for GCs 
from Gen A, the Milky Way, M31, and the Large Magellanic 
Cloud ( LMC). For th e Cen A data, aperture corrections as 
in Hilker et al.l (12007!) have been comp uted for the cluster 
masses and M/Lv (|Mieske et al.ir2008D . The M/Ly values 



^ These isochrones are based onlBertelli et al.1 (|1994| ). but use 
the AGB treatment as in iGirardi et al.1 (|2000D . 



for mass-independent cluster models (the 'canonical' mod- 
els from KL08) are overplotted, and fail to reproduce a 
large part of the data. When considering the metallicities 
and errors of the data (not shown in Fig. [J), only 39% of 
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log Mass [M,r 



Fig. 1. Cluster evolution in the {Af, M/Ly}-plane for mctallicities Z = {0.0004,0.004} ([Fe/H] = {-1.7,-0.7}) 
and dissolution timescales to = {0.3,1,3,10,30} Myr. Solid lines represent evolutionary tracks for initial masses in 
the range = 10^ — 10* Mq with 0.5-dex intervals, while the dotted curves denote cluster isochrones at ages 

t = {4, 8, 12, 19} Gyr. Dots denote the onset of the preferential loss of low-mass stars for each evolutionary track. 
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Fig. 2. Mass-to-light ratio versus mass for globu- 
lar c lusters in di ff erent galaxies. The Cen A data is 
from iMieske et al] ([2008f ). all other data is taken from 



iReikuba et al.l ( 20071 ). The canonical (mass-independent) 
models are overplotted as solid, dashed and dotted lines for 
Z = {0.0004,0.004,0.008}, respectively. 



the observed GCs can be covered within their 1-cr errors if 
the canonical models are use43. 

Although the data show quite some scatter, there are 
indications for a trend of increasing M/Ly with mass, such 
as the lack of low-Af / L clusters for log AI > 6.3. 1 argue that 
this is the same behaviour as can be observed in the cluster 
isochrones including the preferential loss of low-mass stars 
in Fig. [Tl implying that the increase of M / Ly with mass 
corresponds to a decreasing effect of energy equipartition. 

3.3. Explaining tlie mass-to-light ratios of globular clusters 

Figure [3] compares the model cluster isochrones (t — 
12 Gyr) to the GC data for Cen A, the Milky Way, M31 



^ SSP models (e.g., iBruzual fc CharlotI l200g : 
lAnders fc Fritze-v. AlvenslebenI 120031 ) all predict M/L = 2— 
4 Mq for GC metaUicities, a Kroupa IMF and t = 12 Gyr, 
comparable to the 'canonical' models used here. 



to range (Myr) 


Galaxy 


Z = 0.0004 


Z = 0.004 


Z = 0.008 


Cen A 


< 5—5 


< 1—2 


0.2 2 


M31 


< 1—10 


< 0.5—2 




MW 


1—20 


0.7—8 


0.6—0.6 


LMC 


< 3—20 







Table 1. Required dissolution timescale ranges for glob- 
ular clusters of three metaUicities in four studied galaxies, 
as derived from the cluster samples. Due to possible in- 
completeness at low masses and high AI/L (see Fig. [3|), 
all upper limits represent minimum values and lower lim- 
its often represent maxima. Limits that do not suffer from 
incompleteness are shown in boldface. 



and the LMC. The different colours represent three mctal- 
licities, and coloured model lines belong to data points of 
the same colour. 

For any galaxy and metallicity, the data cover an area in 
the {Af, Af/Ly {-plane that can be spanned by two model 
curves of different dissolution timescales. Left curves de- 
note upper limits, while right curves represent lower limits 
for the dissolution timescale ranges in which clusters are 
observed. These limits are chosen as such that they encom- 
pass the data points. Contrary to the sparse coverage of 
the data by the canonical models (see Fig. [5]), it is shown 
in Fig. [3] that 92% of the data can be explained within 
their l-tr errors by using the new models (KL08) that ac- 
count for dynamical effects. The remaining 8% has too high 
M/L to be explained by stellar population models, unless 
their observed ages or metaUicities are underestimated. 

The minimum and maximum dissolution timescales that 
are required to explain the observations are summarised 
in Table [TJ All value s fall within the phy sically reasonable 
rang e of 10^—10® yr (jl amers et al.ll2005bl ) . For each galaxy, 
a broad range of dissolution timescales is required to ex- 
plain the observed range of M/L. This is not surprising, 
since the observed clusters are located at various galacto- 
centric radii and thus experience different tidal dissolution 
strengths. Regardless of this spread, there is a clear trend of 
decreasing required dissolution timescale with metallicity. 
This is likely to be related to the radial metallicity gr adient 
observed in galaxies (first established bv lSearljll97lh . with 
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metal-rich clusters at small galactocentric radii and thus at 
short dissolution timescales. Another trend is that of de- 
creasing required dissolution timescale with galaxy mass. 
Again, this is not surprising, as more massive galaxies can 
have stronger tidal fields and thus give rise to more rapid 
cluster dissolution. 

4. Discussion 

In this Letter, I have shown that the hitherto unexplained 
discrepancy between observations and models of the mass- 
to-light ratios of globular clusters can be explained by dy- 
namical effects. The preferential loss of low-mass stars due 
to energy equipartition gives rise to M/L evolution that 
depends on cluster mass, contrary to what is assumed in 
canonical cluster models. This is confirmed by the appli- 
cation of models that include dynamical effects to the GC 
populations of Cen A, the Milky Way, M31 and the LMC. 

Without the preferential loss of low-mass stars, current 
stellar population models cannot explain mass-to-light ra- 
tios below 2 M© Lq^ for metallicity Z = 0.0004 and below 
2.8 Mq Lq^ for Z = 0.004. As becomes clear from Fig. |3l 
this would leave half of the cluster sample in Cen A and 
most of the Milky Way sample unexplained. Accounting for 
the effects of energy equipartition increases the explained 
percentage of the observations from 39% to 92%. 

The dissolution timescales required to explain the ob- 
served GC samples lie within the physically reasonable 
range of = 10^ — 10^ yr. Observed trends of decreasing 
dissolution timescale with galaxy mass and metallicity are 
as expected when considering the strength of tidal dissolu- 
tion and the radial metallicity gradient in galaxies. 

The dependence of M/L on cluster mass (and thus on 
luminosity) implies that photometrically derived masses 
using canonical models may be strongly overestimated 
(KL08). The results presented here underline the impor- 
tance of accounting for dynamical effects when modeling 
clusters or interpreting observations of (globular) clusters. 
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Fig. 3. Comparison of t = 12 Gyr model cluster isochrones in the {M, M/iy}-plane to globular clusters in four 
different galaxies. Clusters with metallicitics Z = {0.0004, 0.004, 0.008} are denoted with blue diamonds, green triangles, 
and red squares, respectively. Model curves for these metallicities are shown in the same colours, with dotted, dashed 
and solid lines, respectively. In all cases an age of t = 12 Gyr is assumed. For each galaxy, a dash-dotted line of constant 
luminosity represents the fai ntest cluster that is covered by the models , illustrating that the samples are magnitude- 
limited. Metallicities are from lBeaslev et al.l (I2008L Cen Al . iHarrid (fl99l Milky Wavl. lDubath fc Grillmai"^ ([TQQTI . M31) 
and iMackev fc Gilmord (|2004l LMC). To prevent crowding, error bars are only shown for clusters that fall outside the 
range covered by the models. 



